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We consider quasinormal modes with complex energies from the point of view of the theory of 
quasi-exactly solvable (QES) models. We demonstrate that it is possible to find new potentials 
which admit exactly solvable or QES quasinormal modes by suitable complexification of parameters 
defining the QES potentials. Particularly, we obtain one QES and four exactly solvable potentials 
out of the five one- dimensional QES systems based on the sl{2) algebra. 
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Introduction.— Quasinormal modes (QNM) arise as perturbations of stellar or black hole spacetimes They 
. . . are solutions of the perturbation equations that are outgoing to spatial infinity and the event horizon. Generally, 
these conditions lead to a set of discrete complex eigenfrequencies, with the real part representing the actual frequency 
3S of oscillation and the imaginary part representing the damping. QNM carry information of black holes and neutron 
stars, and thus are of importance to gravitational-wave astronomy. In fact, these oscillations, produced mainly during 
^ , the formation phase of the compact stellar objects, can be strong enough to be detected by several large gravitational 
Q ' wave detectors under construction. Recently, QNM of particles with different spins in black hole spacetimes have also 
received much attention. 

Owing to the intrinsic complexity in solving the perturbation equations in general relativity with the appropriate 
boundary conditions, one has to resort to various approximation methods, eg., the WKB method, the phase-integral 
method etc., in obtaining QNM solutions. It is therefore helpful that one can get some insights from exact solutions 
I— I in simple models, such as the inverted harmonic oscillator and the Poschl- Teller potential. Unfortunately, the number 
: of exactly solvable models is rather limited. 
I ' Recently, in non-relativistic quantum mechanics a new class of potentials which are intermediate to exactly solvable 
ones and non-solvable ones has been found. These are called quasi-exactly solvable (QES) problems for which it is 
^ possible to determine analytically a part of the spectrum but not the whole spectrum [^i iSi, i4i, iSj . The discovery of 
this class of spectral problems has greatly enlarged the number of physical systems which we can study analytically. 
In the last few year, QES theory has also been extended to the Pauli and Dirac equations. 

In we have considered solutions of QNM based on the Lie-algebraic approach of QES theory. There we demon- 
strate that, by suitable complexification of some parameters of the generators of the sl{2) algebra while keeping the 
Hamiltonian Hermitian, we can indeed obtain potentials admitting exact or quasi-exact QNMs. Such consideration 
has not been attempted before in studies of QES theory. Our work represents a direct opposite of the work in 01, 
where QES real energies were obtained from a non- Hermitian PT-symmetric Hamiltonian. 
\ QES Theory.— Let us briefly review the essence of the Lie-algebraic approach 0, Q to QES models. Consider 
■ a Schrodinger equation Hijj = Eip with Hamiltonian H = — + V{x) (dx = d/dx) and wave function ^lJ{x). Here 
— . ' X belongs either to the interval {—00,00) or [0,oo). Now suppose we make an "imaginary gauge transformation" 
, on the function -0: ^'(2^) = xi^)^~^^''\ where g{x) is called the gauge function. For physical systems which we are 
^ ' interested in, the phase factor exp(— 5(2;)) is responsible for the asymptotic behaviors of the wave function so as to 
, ensure normalizability. The function xi^) satisfies a Schrodinger equation with a gauge transformed Hamiltonian 
' Hg — e^He^^. Suppose Hg can be written as a quadratic combination of the generators J° of some Lie algebra 
5h ] with a finite dimensional representation. Within this finite dimensional Hilbert space the Hamiltonian Hg can be 
. - - 1 diagonalized, and therefore a finite number of eigenstates are solvable. Then the system described by H is QES. For 
one-dimensional QES systems the most general Lie algebra is sZ(2), and Hg can be expressed as 

Hg=Y^ CabJ^J^ + + ^'^^^ constant , (1) 

where Cab, Ca are taken to be real constants in [1, [1]. The generators J° of the sl{2) Lie algebra take the differential 
forms: — z^dz — nz , J° = zdz — n/2 , J~ = dz {n = 0, 1, 2, . . .). The variables x and z are related by some 
function to be described later, n is the degree of the eigenfunctions x, which are polynomials in a (n+ l)-dimensional 
Hilbert space with the basis (1, z, z^, . . . , z"). 

Substituting the differential forms of J" into Eq. ([1]), one sees that every QES operator Hg can be written in 
the canonical form : Hg = —Pi{z)d^^ + P'^{z)dz -\- ^2(2), where Pfc(z) are fc-th degree polynomial in z with real 
coefficients related to the constants Cab and Ca- The relation between Hg and the standard Schrodinger operator H 
fixes the required form of the gauge function g and the transformation between the variable x and z. Particularly, 
X = dy/ \JPa{v)- Analysis on the inequivalent forms of real quartic polynomials P4 thus give a classification of all 
s;(2)-based QES Hamiltonians If one imposes the requirement of non-periodic potentials, then there are only 

five inequivalent classes, which are called case 1 to 5 in Q. 
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Our main observation is this. If some of the coefficients in Pk{z) are allowed to be complex while keeping V{x) 
real, then all the five cases classified in Q can indeed support QES/exact quasinormal modes. We shall illustrate this 
using one of the cases below. 

QES QNM.- We consider Case 3 in which corresponds to Class I potential in Turbiner's scheme [3]. There 
are two subclasses in this case, namely. Case (3a) and (3b). We shall present the analysis of QNM potential for case 
(3a). The other case turns out to give the same potential with a suitable choice of the parameters. The potential in 
case (3a) has the form (in this paper we adopt the notation of Q): 



V{x) = Ae^^"" + Be^"" + Ce"^^ + De'^^^ 



(2) 



where x £ (— oo, oo) and v is a positive scale factor. Note that V{x) is defined up to a real constant, which we omit 
for simplicity, as it merely shifts the real part of the energy. This remark also applies to the other cases. V{x) in 
Eq. ^ reduces to the exactly solvable Morse potentials when A = B = 0, oi C = D = 0. This potential is QES when 
the coefficients are related by 
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c = i^i!i±l)^ I 

0,1,2... 
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Here 6, c, d are arbitrary real constants. For each integer n > 0, there are n 
the in + l)-dimensional QES subspace: 



1 exactly solvable eigenfunctions in 



^n{x) = exp 



2v 
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Here Xn{z) is a polynomial of degree n 'm z — exp(y^x). To guarantee normalizability of the eigenfunctions, the real 
constants 6, c, d, v and n must satisfy certain relations [![. 

We want to see if we can get QNM solutions if we allow some parameters to be complex, while still keeping the 
potential V{x) real. This latter requirement restricts the possible values of the parameters, and hence the forms of 
QES potential admitting quasinormal modes. For the case at hand, we find that one possible choice of values of 6, c 
and d is: 



b = ib, c = — (n + l)i^, d ^ d, b,d : real constants 



The potential Eq. ^ becomes 

Vn{x) = 

and the wave function Eq. (j4]) becomes 
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Vix) approaches ^oo as a; ^ ±oo respectively: it is unbounded from below on the right. For small positive d and 
sufficiently large n, V{x) can have a local minimum and a local maximum. In this case, the well gets shallower as d 
increases at fixed value of n, or as n decreases at fixed d. We emphasize here that for different value of n, each V(x) 
represents a different QES potential admitting rt+ 1 QES solutions. Since Vix) — )■ oo as a; — > — cx3, the wave function 
must vanish in this limit. This means d > from Eq. ([7|). For the outgoing boundary condition, we must take 5 > 0. 
Before we go on, we note here that there is another possible choice of the parameters, namely, 



'6, c = (n + l)j/, d = id, b,d : real constants 



(8) 



However, this choice leads to a potential related to Eq. ^ by the reflection x i~-> —x. Hence, we will only discuss the 
potential in Eq. ^ here. 
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To see that the wave functions ipn{x) do represent quasinormal modes, we determine the corresponding energy E„. 
This is easily done by solving the eigenvalue problem of the polynomial part of the wave function. From the 

Schrodinger equation we find that for n = 0, the energy is Eq = —u/A — ihd/2v. This clearly shows that the only QES 
solution when n = is a QNM with an energy having a negative imaginary part (recall that h,d,v > Q). For n = 1, 
we have two QES solutions. Their energies are Ei = —hv/A — ihd/2v± \/v^ ~ ibd. Again we have two QNM modes. 
One can proceed accordingly to obtain n + 1 QNM modes with higher values of n. However, for large n, computation 
becomes tedious, and one has to resort to numerical means. 

One can do the same for the other four cases listed in [5j. It is found that these cases admits exact QNM solutions. 

Summary.— To summarize, we have demonstrated that it is possible to extend the usual QES theory to accom- 
modate QNM solutions, by complexifying certain parameters defining the QES potentials. We found that the five 
s^(2)-based QES systems listed in Q can be so extended. While one of these cases admits QES QNM, the other four 
cases give exact QNM solutions. It is hoped that our work would motivate the search of many more exact /quasi-exact 
systems of QNM in QES theories based on higher Lie algebras, and in higher dimensions. 
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